Multi-particle composites in density-imbalanced quantum fluids 
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We consider two-component one- dimensional quantum gases with density imbalance. While genet- 
ically such fluids are two-component Luttinger liquids, we show that if the ratio of the densities is 
a rational number, p/q, and mass asymmetry between components is sufficiently strong, one of the 
two eigenmodes acquires a gap. The gapped phase corresponds to (algebraic) ordering of (p + q)- 
particle composites. In particular, for attractive mixtures, this implies that the superconducting 
correlations are destroyed. We illustrate our predictions by numerical simulations of the fermionic 
Hubbard model with hopping asymmetry. 



Thanks to recent advances in experimental techniques 
of dealing with cold gases, it is now feasible to engi- 
neer one-dimensional (ID) quantum fluids by confin- 
ing atoms in cigar-shaped traps with tight radial con- 
finement By devising an appropriate optical lat- 
tice it is also possible to construct a weakly coupled 
array of such ID "tubes," thus allowing one to study 
the dimensional crossover from ID to three dimensions. 
A number of ongoing and planned experiments deals 
with two-component mixtures atoms of either statistics, 
i.e. Fermi-Fermi (FF), Bose-Bose (BB) or Bose-Fermi 
(BF) mixtures |2] ■ Most of recent theoretical work dealt 
with equal-density mixtures, where a rich phase diagram 
containing both gapped and gapless phases was found 
0) 0j 01 ■ I^or mixtures with unequal densities the ground 
state is generally found to be a two-component Lut- 
tinger liquid [1, ll, 0]- For attractive FF mixtures, su- 
perconducting correlations dominate, thus making the 
ground state a ID analog of the long-elusive Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phase, as confirmed 
both by Bethe Ansatz calculations for integrable models 
and numerical simulations The case of unequal 
mass mixtures — where integrable microscopic models are 
not available — has been studied analytically by means of 
effective field theory 0, 0], and numerically by Monte 
Carlo 111 and time-evolving block decimation (TEBD) 
[lo| methods. A common result which emerges is that 
for strong enough mass asymmetry and/or strong enough 
attraction the system collapses, while for moderate mass 
asymmetry and non-zero density imbalance the ground 
state is again a gapless two-component Luttinger liquid 
with an FFLO-type algebraic order. 

In this Letter we study a generic two-component ID 
mixture with density imbalance within the harmonic fluid 
approach ( "bosonization" ) . We reveal a generic mecha- 
nism which, for a certain relation between the densities, 
opens a gap in the excitation spectrum and completely 
destroys superconducting correlations. We concentrate 
on the properties of the FF mixtures, but our predic- 
tions are applicable to BF and BB mixtures with minor 
modifications. Our findings might also be relevant to spin 
ladder materials with non-equivalent chains in high mag- 
netic fields. We further corroborate our predictions by 



DMRG simulations 12] of a Hubbard model with hop- 
ping asymmetry. 

Consider the mixture of two sorts of fermionic atoms, 
which we label by a pseudo-spin index a =t)i- In the 
bosonization approach we introduce for each species a 
pair of scalar fields <f>a{x) and 6a{x) which vary slowl y on 
the scale of n^'^ , where rto- are the average densities [13|. 
Using the Haldane construction we write for the field op- 



erators '^\{x) 



where the summation over s runs over odd integers s. 
14| For the density operator, no-, this leads to : 
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One of the advantages of the Haldane representation 
^ is that an effective low-energy Hamiltonian can be 
written solely in terms of (jja- and Ho- [l3| . In the non- 
interacting case it is given by Hhcc ~ ^o(0t) + '^o{4'l)i 
where : 
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where v^r are Fermi velocities and K^^ = I the so-called 
Luttinger parameters equal to one in the free case. In 
presence of density-density interactions, Jdxdx' Uaa'{x — 
x')na{x)n'^{x'), Eq. ([2]) is modified in several ways. First 
of all, the s = terms of Eq. ([1]) give rise to an acoustic 
coupling : 



Hi= g dx (aa;0|)(9^0|) 
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where 5 is a forward scattering constant for the interspin 
interactions. More importantly, higher harmonics of Eq. 
^ generate the terms of the form : 
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Here Gs,s' and Gs,s' are (non-universal) amplitudes, and 
hp — TTUa- are Fermi momenta. Since the separation of 
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fast and slow variables is inherent in the bosonization 
treatment, one has to discard in (|4|) the terms which os- 
cillate on the lengthscale ^ kp^ . Strictly speaking, Eq. 
^ is only perturbative in g. On the opposite, Eq. ^ is 
assumed to retain its functional form even in presence of 
generic same-spin density-density interaction, with both 
velocities and Luttinger liquid parameters renormalized 
by interaction terms beyond Eq. Q and various irrel- 
evant operators, e.g. band curvature On a phe- 
nomenological level, we can assume Eq. ([2]) (where, in 
general, 7^ 1) as coming from an underlying micro- 
scopic model, with Eqs. ([3]) and (|4|) regarded as pertur- 
bations. 

Equation (jj]) suggests considering generalized coni- 
mensurabilities of the form 



(5) 



where p and q are relatively prime integers. Notice that 
this condition does not imply the presence of a lattice : 
we only require the densities to be commensurate with 
each other. Eq. ([5]) selects from Eq. ^ the terms with 
s/s' =p/q, and the Hamiltonian ([4]) reduces to 



(6) 



H2 ~ G dx cos 2 {p<p^{x) — q4>i{x)) , 



where we only keep the lowest order term, since the scal- 
es 



24| 



ing dimension of the operator cos s(p is s' 

We now assume that the densities are commensurate 
via jSJ, and analyze the model H. = Ho(0t) -I- Hoi4>i) + 
Hi + H2 defined by ©-([S]) and Since in general this 
model is not exactly solvable, the nature of the phases 
can, in principle, be determined by an approximate renor- 
malization group (RG) procedure. Rescaling the fields 
via (/)| = p4>^ and 0| = q(j)\ the model is brought to 
the form considered in Ref. [5[, where an RG procedure 
has been carried out including the renormalization of ve- 
locities Va [see also Ref. in the fermionic language]. 
For large velocity asymmetry and strong attractive (re- 
pulsive) interactions the system was always found to col- 
lapse (phase separate). Barring such an instability, two 
regimes were found, corresponding to the cosine operator 
([S]) being relevant or irrelevant in the RG sense. 

In the regime where the cosine operator Eq. © is ir- 
relevant we are left with a bilinear Hamiltonian ©-([S]), 
which is diagonalized by appropriate linear combinations 
of the fields ^161]. As a result one obtains an effective 
theory Ha which features two decoupled massless fields 
If 1^2 with corresponding velocities vi^2 and Luttinger pa- 
rameters Ki^2- Ha = Hoifi) + Ho{f2)- For an attrac- 
tive FF mixture, such a theory describes a ID analog 
of the FFLO phase : all correlations are algebraic in real 
space and the pair correlation function oscillates with the 
FFLO momentum Qfflo = \kp — kp\. 

Another regime corresponds to the case where the co- 
sine in Eq. ^ is relevant in the RG sense. Then the 



system has a massive mode (pa and a massless mode, (j)b- 
The effective theory, Hb, can be written as : 



HB=HsGi<l)a)+HQi^b) , 



(7) 



where HsG{4>a) — Hoifj^a) + G Jdx cos2-\/20a is the sine- 
Gordon model for the field (pa- The Eq. ([7]) is charac- 
terized by two mode velocities Va.t and two Luttinger 
exponents Ka.t, with Ka ^ 1, so that (pa is pinned by 
the minimum of the cosine operator in ([7]) . Closed- form 
expressions for the parameters of Eq. ([7]) can be easily 
obtained in several limiting cases. Indeed, for — 
the exact eigenmodes oiH = Ho{(j)]) +Ho{(pi) + H2 are : 

(pa = {p(p^ - q(Pi)/V2 , 
^b^{qK^^^+pK^^^)/V2, (8) 

with the Luttinger exponents 

Ka = {p^K^ + q^K{) /2 , KaK^K^ . (9) 

Notice that for higher-order commensurabilities (larger p 
and q) smaller values of K are required for (pa to acquire 
a gap, cf. Eq. (ID). 

Deep in the massive phase one can make a crude ap- 
proximation to the cosine operator in H2 by replacing it 
with a mass term cx (p0| — q(pi)'^ ■ This leads to: 
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Kb = -K^Ki^fUyUl 



{p^K^+q^K^f 
p^K^vi +q'^Kiv^ 



(10) 
(11) 



which reduces to (O for = Wj- 

We now turn our attention to an interpretation of the 
theory ([7]), focusing on the novel regime with p > 1. 
Obviously, excitations corresponding to eigenmodes (pa,b 
carry both spin and charge. Furthermore, these excita- 
tion correspond to multiparticle states in terms of the 
original t, i particles since a particle of the species a cor- 
responds to a 277 kink of the field (pa [l^ • To gain further 
insight to the structure of the massive phase we consider 
its correlation properties. We classify operators 0{x) 
according to whether the asymptotic decay of the corre- 
lation functions (0(0)0^ (x)) for a; — > cxD is exponential, 
(X e""^"!^!, or algebraic, oc jxl"^"*^. For equal densities, 
p = q — 1, the dominant algebraic order [i.e., the small- 
est decay exponent ao ) is found among the two-point op- 
erators : the superconducting fluctuations, Og 



and charge density wave, Ocdw = Y^a.a' 'WRa^'jci'La' , 
and spin density wave, 0%^^ = Y.a,a' ^^^'^R.j^^aa'^La' ■ 
Here ipL,Ra are left- and right-moving fermions, respec- 
tively, and a" are the Pauli matrices. 

The case p ^ q \s markedly different : First of all, 
the superconducting correlations described by Os always 
decay exponentially, and so do the x- and y-components 
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of OgDw- ^'^^ CDW and SDW^ operators we write 
OcDW =0[R + OiR andOi 
the auxilliary operators 
find for a; — s- cx) 



C'lr)' where 



i'Lcr^R'y- Using ^ we 



where = q^Kbl2Kl, 6| = pKi^/Ka and = 

i\/20a6T^ Likewise, for the J,-species the exponent 

is Qfj = p^Kb/2Kl and the amplitude is A{bi) witli 
fej = qKi/Ka- The ampHtudes A(6) depend exponen- 
tially on & : logA(&) oc 1/6* [l3|- We thus see that cor- 
relations of OcDW and OgDw both given by a super- 
position of two power-laws (|12p with exponents a|^| — 
where the slower the decay, the smaller (exponentially 
smaller) is the corresponding amplitude. 

Given the massive mode in the form ([5]) with p ^ q, we 
construct a compound operator Op+g — which has 

algebraically decaying correlations. Specializing for the 
lowest order commensurability ([5]) with p = 2 and 9=1, 
this corresponds to a "trimer" operator ©2+1 — ^j^j^j- 
For fermionic t component the corresponding decay ex- 
ponent a2+i = {Kb/2Kl + 2Kl/Kj,) /2. We thus see 
that in this particular case the dominant correlations 
in the massive phase are the 2kp density waves ([12]) 
for < l/-\/3, and the "trimer" correlations Op+g for 
a| > l/\/3- We stress that the competition between Olr 
and Op+q is generic, in a sense that it holds irrespective 
of the statistics of t- and |-particles both on the lattice 
and in the continuum. 

Microscopies. — We now focus on the following ques- 
tion: Is there a microscopic model whose low-energy ef- 
fective theory would be given by Eqs. ([5]) and (O? 

We start from constructing such a model explicitly in 
the weak-coupling regime with respect to the interspecies 
interaction. Namely, for the J, component we take non- 
interacting fermions (or, equivalently, Tonks bosons) of 
the mass and (Hnear) density n;, so that = 1 
and vi — -Kni/mi. For the ] species we take a dipo- 
lar Bose gas which is known to be a Luttinger liquid 
with 7r[6C(3)n|ro]^"'^/^ as u^tq — > 00 18]. Here 

tq = m^cP 12-K is the effective Bohr radius associated 
with the dipole moment d and C, is the Riemann zeta 
function. We thus see that for n^r^ — p''7r-^/6C(3) we 
have = 1/p^. Furthermore, Galilean invariance fixes 
the product W| ~ ttu-^ /m| [14| . Constraining the den- 
sities via ([5]) withp > g = 1, and assuming m| = pm^ we 
have both U| = and p^K^ = q^K^ = 1 by construction. 
Now, coupling the t and | species via, e.g., a short-range 
interaction U J dxnij{x)ni^{x) with infinitesimal U gener- 
ates the terms of the form ^ and ^ with g = U /tt^ . 
The eigenmodes of the system are then given by Eq. ([5]) 

and a direct calculation yields Ka ~ P^K^ {l + 3-^^ 

0(52),andX,=p2^2^^ (l_g:^^)+0(g2)^ ^^^l^^^ 



immediately see that having U < Q yields Ka < 1 and 
hence drives the system to the gapped phase, where the 
gap A is exponentially small: In A ~ — const/?/. A simi- 
lar construction can easily be effected for an FF mixture 
on a lattice. In this case we take for the t component, 
e.g., a model with finite-range interactions [l9j . 

In the example above we engineer the theory ([5]) and 
(l7|) by coupling a majority of light and non-interacting J, 
species to the minority of heavy particles t, which have 
strong repulsions among themselves. Such a construc- 
tion is somewhat ad hoc, and requires fine-tuning. A 
much more natural alternative is provided by a simple 
observation : even purely local interspecies coupling U 
generates long-range effective interactions in higher or- 
ders of perturbation theory. Thus, having finite U and 
m| ^ m| should be sufficient to divert the RG flow to- 
wards the theory ([7]). In this case we expect particles of a 
light minority component to provide an effective coupling 
between heavy particles of majority species. 

To this end we consider an asymmetric attractive {U < 
0) Hubbard model 



aH 



E 



h.c 



U 



E^'T"-a ' (13) 



where Cia- annihilates a fermion with spin a on a site i € 
[1, L] of a chain lattice of length L, (ij) stands for pairs of 
nearest neighbor sites, hi^ = cl^Cia, and ta are hopping 
amplitudes for the spin-up and spin-down components. 
For rj = ti/t^ — 1 the model is solvable by Bethe 
Ansatz techniques even in presence of density imbalance. 
For rij 7^ ni the ground state is of the FFLO type [§], 
which, in present langua ge, corresponds to a gapless fixed 
point theory Ha — see |20| for a detailed discussion. 

For unequal hopping amplitudes, 77 ^ 1, the model 
(|13p is no longer integrable, and we resort to numerical 
simulations using DMRG technique [l2| • We use lattices 
of up to L = 80 sites with open boundary conditions 
and DMRG truncation of up to TVs = 400 states and 
check that (i) the discarded probabilities amount to no 
more than 10~^, and (ii) the results are stable with re- 
spect to Ns- We calculate single-particle density matri- 



ces pa{x) = (cl/2,, 



and pair-pair correlations 

r(a;) — (^L/2'^l/2+a;)j where Vj = Cj^Cji is the lattice 



de- 



version of the superconducting operator Og and { 
notes an expectation value over the ground state. 

Fig. [1] shows typical results for the pair-pair cor- 
relations r(a;). We find that for small enough hop- 
ping asymmetry, 77 > rjci, the long-distance decay of 
both single-particle (not shown) and two-particle corre- 
lations is consistent with the FFLO-type laws r(a:) oc 
cos((3fflo 2^)|a;|~''' and Pa{x) oc cos(fc^a:)|a;|~'^. On the 
contrary, once the hopping asymmetry exceeds some crit- 
ical value and the densities are commensurate via ([5]). 
the power-law decays change to exponentials, namely 
T{x) oc e~l^l'^ cos((5fflo a;)|a;|~''' and likewise for Pa{x), 
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FIG. 1: Superconducting correlation function r(a::) for the 
asymmetric Hubbard model (|13|l for ni = 2n| = 3/10: rj — 1 
(black squares) and rj — 0.3 (red circles). Shown by blue 
triangles is r(a::) for rj — 0.3 and 'incommensurate' densities 
n| = 17/80 and = 29/80. Hubbard coupling is [/ = -5tf 
and the system size L = 80. Lines are guides to the eye. Inset: 
Fourier transform of r(a;), same color coding. Arrows indicate 
the characteristic momenta: kp — kp, kp and kp + kp, respec- 
tively, for the 'commensurate' densities — 2n| = 3/10. We 
stress that in all these simulations the density distributions for 
both components are uniform apart from Friedel oscillations 
induced by the open boundary conditions. 



thus unequivocally signalling the presence of a gap. Vio- 
lating the relation ([5]) destroys the gap, and the correla- 
tion functions decay algebraically again. 

The inset in Fig. [1] shows the superconducting corre- 
lation function in momentum space. Compared to the 
Hubbard limit {rj = 1), we see that the mass imbalance 
[t] — 0.3) leads to an overall broadening of the distribu- 
tion, which now extends well beyond kp + kp. In addi- 
tion, the opening of the gap at commensurate filling de- 
pletes the superconducting correlation at small momen- 
tum. Detailed investigation of the asymmetric Hubbard 
model ([T^ is beyond the scope of this Letter and will be 
reported elsewhere [21 1. 

Conclusions and outlook. — Summarizing, we have re- 
vealed a generic mechanism of opening a gap in two- 
component quantum fluids with density imbalance in one 
spatial dimension. The gapped phase appears once in- 
teractions and mass asymmetry between components is 
strong enough, and the densities satisfy Eq. ([5]). Depend- 
ing on the microscopic details, the system develops quasi- 
long range ordering of either 2kp density waves or of pe- 
culiar (p + (7)-particle composites. The proposed mech- 
anism applies to mixtures of particles of either statis- 
tics, and does not require the presence of a lattice. Ex- 
perimental signatures of the proposed state include (i) 
the disappearance of the superconducting ordering, and 
(ii) appearance of the (p -I- g)-particle composites, which 



can be detected, e.g., by noise correlation measurements 
in the time-of-flight absorption imaging using the tech- 
niques discussed in Ref. [2a |. 
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[24] For lattice models there is an additional possibility that 
pn^ + qni — integer/ao, where ao is the lattice constant. 
In this case an additional term of the form cos(p<jf>| + 



would appear in the Hamiltionian. The effects due to such 
term are similar to those due to Eq. 



